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Abstract 

In this letter we derive new modes for massive Dirac field in the background of a higher¬ 
dimensional Schwarzschild black hole. We use in our approach the Cartesian gauge defined 
in local frames. We work in the context of ADD-like theories, assuming that the Standard 
Model fields (fermions and bosons) live only on a (3 + l)-dimensional brane and the extra 
dimensions of space can be large in size. 
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1 Introduction 

One of the main goals of today’s physics is to develop a unified picture of nature in which all 
four of the known forces (gravitational, electromagnetic, week and strong nuclear force ) will be 
unified. During the last decades many ideas were proposed for the unification of gravity with the 
rest of the forces. Many of these ideas are based on the possible existence of extra dimensions 
of space. Although at the beginning it was thought that all extra dimensions must be very 
tiny in size and compact, the work of Arkani-Hamed, Dimopoulos and Dvali(ADD)[l, 2, 3] 
have showed that these (linear) extra dimensions could have submillimeter and even millimeter 
length scales. In these Large Extra Dimensions scenarios a solution for solving the hierarchy 
problem was proposed[l], thus lowering the fundamental energy scale of gravity (Plank scale) 
up to the order of TeV. 

In the ADD brane-world scenario all the Standard Models particles (fermions and bosons) 
are constrained to live only on a (3 + l)-dimensional brane (our observable 4-D Universe), while 
in the bulk only gravity can propagate and possibly other unknown scalar fields. Thus one 
possible way for detecting the extra dimensions is to study how gravity affects the other fields 
that live on the brane. This can be done, for example, by studying higher dimensional black 
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holes. The process of higher-dimensional black holes creation at TeV scale energies was in the 
past decade the subject of many articles[4, 5, 6, 7, 8, 9]’ [10, 11, 12, 13, 14, 15, 16] . It is assumed 
that once such objects are formed they will quickly evaporate through Hawking radiation[17]. 
This evaporation process of a higher-dimensional black hole was also the subject of many studies 
in the literature[18, 19, 20, 21, 22, 23, 24, 25, 26]. 

In the vast majority of these studies the field equations, that describe the motion of the 
particles in the (3 + 1)11 brane induced gravitational background of the higher dimensional 
black hole, where studied manly for massless particles [18, 19, 20, 21, 22], One of the reasons 
for that is that these studies are using the Newman-Penrose[27] formalism in which the massless 
field equations may be treated more easily (for the massive case in this formalism see Ref. [30]). 
However, in order to treat the massive fields (especially the Dirac field), a more standard 
approach, for dealing with the spin of the field, proves to be more efficient. As a consequence, 
the massive fields in the context of higher-dimensional black holes have been so far less studied 
in the literature. 

In this short letter we continue the work of Ref.[28, 29, 30] by solving the massive Dirac 
equation in the background of a Schwarzschild (4 + n) higher-dimensional black hole, using a 
different gauge technique than the ones used there. We chose to work in the so called Cartesian 
gauge proposed some time ago in Ref. [31, 32], The advantage of this gauge is that the separation 
of the angular variables can be done in a very similar way as in the central problems from flat 
space-times. Thus the angular part of the massive Dirac equation is completely solved in terms 
of the usual 4-dimensional angular spinors [34, 35]. For the remaining radial problem we will use 
an approximation technique for finding analytical solutions in the near horizon zone, respectively 
in the far field zone. 

We would like to underline that the solutions found by us in the near horizon zone form 
a new set of modes that have never been discussed elsewhere in the literature before. These 
modes can be expressed as a combination of Whittaker M and W functions and besides other 
relevant physical constants (like the black hole and fermion masses, angular quantum number 
k, etc) they also contain the dependence on the number n of extra space dimensions. Regarding 
the far field modes we have found that they are a combination of Bessel functions of the first 
(J) and second kind (Y). We have thus recovered in the far field zone similar solutions to those 
reported in ref[28, 29]. To our best knowledge all the solutions to the Dirac equation, found in 
the literature, in the far field zone, for the fermions located on the brane, are independent of the 
number of extra-dimensions. Our solutions are also included in this category. This is no surprise 
because an observer placed in the far field zone will see the black hole as a 4-dimensional object 
due to the compactification of the extra dimensions. 

Starting form the above mentioned new set of modes, one can in principle use the matching 
technique to combine the near horizon ones to those from the far field zone in order to calculate 
analytical expressions for the massive fermion absorbtion probability and most importantly for 
calculating the low-energy Hawking radiation (and luminosity) emitted during the Schwarzschild 
phase[4, 36, 37] in the life of a higher dimensional black hole. 

In the next section (2) we present some key features of the Castesian gauge that we are 
using to study the massive fermion fields. After deriving, in the first part of section 3, the 
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general equation for a massive Dirac field constrained to move on the (3 + 1) brane projected 
background of a (4 + to) higher-dimensional black hole, we proceed with our investigation by 
studying the near horizon (section 3.1) and far-held (section 3.2) behavior of spin 1/2 particles. 
In the last section (4) we present our conclusions derived from this research. 


2 Preliminaries 


In this section we will briefly review, following Ref[31, 32] , how one can write down the Dirac 
equation in central charts using the Cartesian gauge [31] which explicitly points out the central 
symmetry of the manifold. The Dirac equation in curved space-times can be defined in a local 
frame chart as [31, 38] 

— rm/j = 0 , (1) 

taking the explicit form [31, 33] 

- mf’ + ^-J=d„(y/=ge£)'Y a ip - ^{ 7 “, S^}t'l.i/t = 0 , ( 2 ) 


Both natural indices, a ,.., gt, v ,... and local indices a,ft ,... take the same values 0,1, 2,3 
and S^~ stands for the generators of the spinor representation of the SL(2,C) group. The 
line element can be expressed in terms of the tetrads e&, e a and the 1-forms = ivdx u as 
ds 2 = rj & pLo a ujP = g^dx^dx 1 '. Using the tretrads defined in ref. [31] one can show that the 
explicit form of the line element in a central static chart with spherical coordinates (t, r, 8, </>) 
depends on three arbitrary functions u(r), v(r) and w(r) and has the following expression[31] 


ds 2 


w(r) 2 



dr 2 

u(r) 2 


2 

, r 0 (d9 2 + sin 2 ddcj) 2 ) 
v(r) z 


( 3 ) 


One of the main advantages of working with the Cartesian gauge is that the Dirac equation 
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can be reduced to a much simpler form with the help of the transformation if —>• vw~^ip. The 
new equation[31] allows us to separate the spherical variables in the same manner done for 
the central problems in flat spacetimes. The positive frequency solutions of energy E to this 
equation can then be written as 


v(r) 


<f>) = —— ^ [/g |W (r)$+ j . )II (g, (f) + / g|B (r)$ mj .,K(fl, ^)l e 


—iEt 


( 4 ) 


where K are the usual four-component angular spinors [34, 39] and /^ K (?') represents the 
radial wave functions. One can show that these functions satisfy a system of radial equations, 
written in compact form as 


1 m w(r) -u(r)J: + K^\ / /+» 

\ u(r)^+K^?X —mw(r ) / \ fE,S r ) 


= E 


( fE.S r ) \ 


V -W 7 ") / 


( 5 ) 


the radial Hamiltonian H r 
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In what follows we will solve this eigenvalue problem for massive fermions constrained to 
move only on the brane. 

3 The Dirac equation in higher dimensional Black Holes 

In theories with extra dimensions of space one can discuss various problems. Among the most 
important ones is the study of black hole production[41, 36, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 

14, 15, 16] and the proprieties such objects may have [48]. Once produced these mini black 
holes will evaporate through Hawking radiation almost immediately. During this process the 
black hole will pass through several phases[4, 36, 37]: (1) the balding phase during which the 
black hole emits mainly gravitational radiation; (2) the spin-down phase , in which the initial 
angular momentum of the black hole it is lost via Hawking radiation and superradiance; (3) 
the Schwarzschild phase when the black hole becomes spherically-symmetric and it’s mass will 
decrease as the Hawking temperature increases due to the emission of radiation; (4) the Plank 
phase when the mass of the black hole is close to the Planck scale, phase for which we need a 
quantum theory of gravitation. 

In our study we will concentrate on the Schwarzschild phase in the context of ADD-like 
theories[1, 2, 3]. The purpose of this section is to study the brane fermion fields by solving the 
massive Dirac equation in the 4-dimensional induced background of the bulk black hole metric. 
For this we will apply the procedure described in section 2. We would like to point out that, 
up to our knowledge, this is the first time when one tries to study the brane fermions using the 
cartesian gauge[31] for solving the Dirac equation. 

A key feature of the ADD-like theories is the assumption that the Standard Model fields 
(fermions and gauge bosons) live only on a 4-dimensional brane, thus they will propagate in the 
gravitational background described by the 4 D induced metric: 

ds 2 = fir) dt 2 —- — r 2 (d0 2 + sin 2 8 d<f 2 ) , (6) 

fir) 

obtained from the higher-dimensional (4+n)-Schwarzschild black hole which has the line elemental, 
40] 

ds 2 = f{r) dt 2 - - r 2 dPt 2 2+n (7) 

The induced metric (6) is obtained form (7) by setting 6i = 7r/2 (i > 2) in the angular part of 
the metric dQ 2 , n gi ven by 

d,n 2 +n = dd 2 +1 + sin 2 9 n+ i { dd 2 + sin 2 9 n [... + sin 2 0 2 (d0 2 + sin 2 9\ d<f 2 )...\ } (8) 

The function /(r), both for the general metric (7) and the induced one (6), has the following 
expression 

/<o—(?r «» 

where n gives the number of space-like extra dimensions and ru represents the horizon radius, 
obtained by applying the Gauss law in a (4 + n)-dimensional spacetime. Now th has a more 
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complicated expression (than the usual Schwarzschild, ro = 2MbhG, radius) which depends 
on the number of extra dimension n, the mass of the black hole Mbh and the fundamental 
higher-dimensional Plank mass M* (related to the standard Plank mass by: Mp ~ M 2+n R n , 
with R the size of the extra dimensions) given by [41] 


1 (Mbh\ ^ ( 8£Pg)\ "+ 1 

7tM* \ M* / \^n + 2 J 


( 10 ) 


Comparing (3) with (6) we can identify the functions 


u(r) = f(r) = 1 - (y-^j + , w(r) = v(r) = y/f(r) = \jl- (y^j (11) 

which enter into the radial Hamiltonian H r given in (5). As done in Ref[42], (for the standard 
4D Schwarzschild black hole) we introduce the Novikov dimensionless coordinate[43] 


x = 



1 € (0, oo), 


( 12 ) 


with the help of which the exact radial eigenvalue problem H r J- = EE in the space of two- 
component vectors E = (/ + , f ~) T , can be rewritten as follows 


(A B\ 

( f + {x) 

U d) 

V /”(*) / 


(13) 


where we introduced the notations n = , e = rnE and 


A 

C 


v? r 

Id k 

2x dx (x 2 + l)v? ’ 


Id k 

2x dx (x 2 + l)v? 



/ 


(14) 


In the next subsections we will solve equation (13) by searching for analytical solutions in 
the vicinity of the black hole’s horizon, respectively far away from it. 


3.1 Solutions in the near-horizon zone 

For finding the modes in the vicinity of the black hole we first need to approximate the operators 
given in (14) by expanding them for small values of x. This can be done using a Taylor expansion 
for x ~ 0, obtaining the following expressions: 


1 _ i 1 
/ x n + 1 

j_ * r±L 

y/J X V n + 1 


1 n + 2 n(n + 2) 3 5 

x + — +0(I) 


i + 


n + 2 


+ 


n 2 — 4n — 12 


96 


x 4 + 0(x 6 ) 


1 


1 


(x 2 + l)v7 x V n + 1 


1 n — 2 2 n 2 — 28n + 36 4 , 

1 + ——x 2 +-—- x 4 + 0(x 6 ) 

4 96 


(15) 
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Keeping only the relevant terms (i.e. those proportional to x, 1/x, respectively any constant 
terms) and neglecting the ones in 0(x 2 ) the system given in (13) can be rewritten as 




y/n + 1 n + 1 \x 


1 n + 2 
- H- ^—x 


fnh ( x ) - ( - 


2dx y/nTl 1 Jnh 


fnh ( x ) = 0 


+ 

ld + ^=|/,i(x)- 


2 dx \J ti +"T 


^ + e 


\/ra + 1 n + 1 \x 


1 71 + 2 

H--—x 


(16) 


fnh ( X ) = 0 


x> 


As it is, the above pair (16) of differential equations can not be solved analytically. However, if 
we perform the transformation [42] 




(17) 


/ 1 -i\ 

\~i 1 / 

we obtain a new pair of differential equations for the new radial wave functions (f^ h - f~ h ) : 

[(A ± C) + i(B ± V)\ /± (x) + [(£ ± P) + i(A ± C)\ ft h (x) = 0 (18) 

The analytical solutions of (18) can be found in terms of Whittaker M and W functions[44] as 
1 ^ .,„._2, 1 


fnh ( X ) = ~l= Cl M a,f}(il x ' Z ) + -^=C 2 W Qt p{i^X 2 ) 
v X y/X 

fnh( X ) = (^P+~1 


'nh\~s \^ 2 J 
1 iy/n + 1 


1 - - - * • ^ 1 


-^=CiM ai/3 (ryx 2 ) + -j=C 2 W a> p(i'yx 2 ) 

_y/X y/X 


X K — Ifl 

where we have introduced the following notations 


x 

7, + a + P ) A=C 1 M a+1 ^(i'yx 2 ) - ]=C 2 W a+1 ^(i-fx 2 ) 

^ J v x V X 


(19) 


1 

a = - — i 


H 2 + K 2 


n + 1 e(n + 2) 




7 = 


n + 1 


(n + 2) 


( 20 ) 


P = 


y(n + 1) + e(n + 2) 
y/n + 1(k — i[i) 


Q = 


e + z(n + l)(| + f) 
Vw + 1 (k — in) 


while Ci and C 2 are two arbitrary constants. These new set of modes are one of the main 
results of our research. 


3.2 Solutions in the far-field zone 

Let us now turn our attention to the far-held zone for which r rjj. This can be done by taking 
the limit x —> oo in the relations (14). By doing so we have found that there is a difference 
between the case n = 0 (no extra dimensions) studied in Ref. [42] and the case when one or more 
extra dimensions are present (n > 1). Practically we have found that the Dirac equation has 
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different particular solutions for the two cases. The solutions found in Ref. [42] for n = 0 are a 
combination of Whittaker M and W functions, while our solutions for n > 1 can be expressed 
in terms of Bessel functions. 

In order to obtain the system of differential equations resulted from (13) for the far-field 
zone, we first need to calculate the following expansions in the asymptotic region of the black 
hole: 


1 1 


1 


7 ~ x 

X ' x 2n+l 

1 _ 1 

rn 1 

1 

V? ~ * 

~ ' 2x 2n+1 

i 


1 

'1 

(x 2 + I)v7 X 

X 


+ 0( l/x 2n+5 ) 


n + 1 

rj.2n+3 

+ 0( l/x 2n+3 ) 

. J_ + 0(1/I 5 ) 


( 21 ) 


which are in fact the Taylor series with respect to 1/x of the operators given in (14). Keeping 
only the terms proportional to x, respectively 1/x and any constant terms, neglecting thus the 
terms of the order 0( 1/x 2 ), we arrive at the following pair of equations for the radial wave 
functions (/+,/“): 


fa(x) - x(e + fi)f a (x) = 0 


which is a system of two Bessel differential equations[44] with the general solutions (in the 
asymptotic zone) given by 


fa ( x ) = y/e + n axJ K+ i (vx 2 ) + bxY K+ i {vx 2 ) 
fa (x) = \Je - n axJ K _i(isx 2 ) + bxY K _i(isx 2 ) 


(23) 


where we denoted by v = ^/e 2 — fi 2 , while a and b are two arbitrary coefficients. 

Because our modes are a combination of Bessel functions of the first and second kind, one 
can not talk about the existence of discrete energy levels like in the case of black holes with no 
extra dimensions[42, 47]. This is again a difference between the classical Schwarzschild black 
holes and their higher-dimensional counterparts. We think that a full investigation of the exact 
radial problem (13) using numerical methods[45, 46, 47] for solving the Dirac equation (in the 
Cartesian gauge) may reveal the existence of discrete energy levels for fermions moving in the 
background of a higher-dimensional black hole. 


4 Conclusions 

In this study we investigated the brane behaviour of massive fermions in the vicinity and far 
away from a (4 + n)-dimensional black hole. This was done by finding new analytical solutions 
to the Dirac equation in the near horizon, respectively in the far-held zone of the black hole. 
We have thus proved once more the advantages of using the Cartesian gauge for investigating 
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the (massive) Dirac equation in different curved backgrounds. 

This research can be continued in a number of ways. One of these is the full investigation of 
the radial problem with the help of numerical techniques in order to find a complete (numerical) 
solution. It would be also interesting to derive the Hawking radiation of massive fermions using 
our new modes given here, and to see to what extent this radiation will differ from the one 
obtained in other studies found in the literature. 
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